CU-TP-951 



QCD hadron spectrum with domain wall fermions 

Lingling Wu a * 

a Department of Physics, Columbia University, New York, NY 10027, USA 



We present the QCD hadron spectrum for the cases of both quenched and two-flavor dynamical domain wall 
fermions. We compare the results obtained using the Wilson gauge action and a renormalization group improved 
gauge action. Finite volume effects and the dependence on the finite extent of the fifth dimension are discussed. 
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1. INTRODUCTION 

The domain wall fermions (DWF) formal- 
ism JlJ uses an extra space-time dimension to 
separate the chiral limit from the continuum 
limit. In this paper, we report a study of the 
hadron spectrum obtained from both quenched 
and dynamical DWF. 

Our conventions are: my is the 4-d bare quark 
mass that explicitly mixes the two chiralities on 
the domain walls, L s is the extent of lattice in 
the fifth dimension, and too is the 5-d bare quark 
mass, which is often called domain wall height. 
The masses in physical units in this paper are 
obtained by using the rho mass to set the scale. 

2. QUENCHED QCD SPECTRUM 

Last year, we reported that with domain wall 
fermions at (3 — 5.7 on an 8 3 x 32 lattice with 
to = 1.65 and L s — 48, toat/to p = 1.42(10) as 
to,/ — » o||. Considering the moderate size of the 
lattice, this value is favorable compared with Wil- 
son and staggered results. However, an L s study 
for 8 3 x 32 lattices at this (3 shows that (to/ — > 
0) fits well to the form of Aexp(—aL s ) + B 
with a non-zero value of B w 0.048, which gives 
m v /a(mf -> 0) « 213MeV at infinite L.g. In 
Figure Fy we show the pion mass data on which 



the conclusion was based. 
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• Wilson m =1.65 8 J x32 

♦ Wilson m =1.65 16 3 x32 
> Iwasaki m =1 .40 8 3 x32 
□ Iwasaki m =1.65 8 3 x32 
A Iwasaki m =1.90 8 3 x32 
< Iwasaki m„=1.65 1 6 3 x32 
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Figure 1. Quenched m^.(m/ — > 0) versus L s for 
the Wilson gauge action (J3 = 5.7, closed sym- 
bols) and the Iwasaki gauge action (j3 = 2.2827, 
open symbols). The solid line is 0.048(9) + 
0.37(5) exp(-0.09(2)L s ), with y?/dof = 2.8/2. 



In order to investigate how much of this 
213MeV pion mass comes from the effects of finite 
volume, we have studied a 16 3 x 32 lattice at (3 = 
5.7, L s = 24. We find m\(m s -> 0) = 0.077(2), 
only 0.021 smaller than the 0.098(7) value we ob- 
tained for 8 3 x 32 at L s — 24. This finite volume 
shift is less than half of the 0.048 infinite L s limit 
for 8 3 x 32 suggesting this non-zero limit of to^ is 
not a finite volume effect. 

To examine this question more carefully, we as- 



2 



sume the effects of L s can be represented by a 
residual quark mass, m les (L s ), and express mi 
to first order in chiral symmetry breaking as: 
ml(V 7 m f ,L s ) = co(V)+cx(V)*(mf+m lea (L s )). 
We find ci(8 3 ) = 4.54(9), ci(16 3 ) = 4.75(3). If we 
require that Co vanishes as V — > oo and assume 
that 16 3 is sufficiently large that co(16 3 ) « 
then either m rcs does not vanish as L s — + oo 
or m Iea (L s ) has a strong, unphysical dependence 
on the spatial volume with m res (24) increasing 
from 0.011(2) to 0.0162(5) as the volume is in- 
creased from 8 3 to 16 3 . Thus, our results appear 
to require one of three unexpected possibilities: i) 
quenched m^ does not vanish in the chiral limit; 
ii) m rcs does not vanish as L s — > oo or iii) m res 
has a strong volume dependence for L s = 24. 

To support our weak matrix element study 
we have also measured hadron masses using DWF 
and the Wilson gauge action at a weaker coupling, 
(3 = 6.0. On a 16 3 x 32 lattice with m = 1.80, 
L s = 16 and nif ranging from 0.01 to 0.04, we 
have(Table 0) m N /m p = 1.37(5), m£ = 0.014(2) 
as mf — > 0, which gives m 7r /a(m/ — ► 0) = 
230(15)MeV. 

From the above discussion, we can see that to 
decrease the pion mass while using DWF and the 
Wilson gauge action, large L s is needed. How- 
ever, the simulation difficulty is proportional to 
L s . It might be helpful to use the renormal- 
ization group improved gauge action of Iwasaki 
which smooths out the gauge fields on the lat- 
ticef). 

We have simulated using the Iwasaki gauge at 
(3 = 2.2827 and c\ = —0.331 which is equivalent 
to quenched (3 = 5.7 with the Wilson action. The 
resulting pion masses are promising. Compared 
with the large m\ (m / — > 0) values using the Wil- 
son action, much smaller values are obtained from 
lattice size 8 3 x 32 as shown in Figure |l|. Studies 
with L s ranging from 16 to 32 and mo ranging 
from 1.40 to 1.90 show that m^ has little depen- 
dence on L s and m^. A simulation at larger vol- 
ume (16 3 x 32) also shows that the finite volume 
effect is small(Figure|l|). These suggest that using 
the Iwasaki gauge action may enable us to study 
quenched DWF at smaller Ls, but we may still 
have a non-zero pion mass at infinite L s . 



Table 1 

Valence extrapolations (a + bm^ 1 ^) for mi, m p 
and mjv- Quenched fits are obtained on a 16 x 
32 lattice with mo = 1.80. Dynamical fits are 
obtained on an 8 3 x 32 lattice with mo = 1.90. 



mass f3 


L s 


a 


b 


quenched: 


m^ o.u 
m p 6.0 

mjv 6.0 


16 
16 
16 


0.014(2) 
0.400(8) 
0.55(2) 


z.yo^ / ) 
2.80(13) 
4.7(3) 


dynamical with 


m ( f dyn) = 0.02: 




mi 5.325 
m p 5.325 
m^ 5.325 


24 
24 
24 


0.325(4) 
1.13(2) 
1.62(3) 


5.23(3) 
1.80(6) 
3.48(13) 


m'i 1.90 
m p 1.90 
m N 1.90 


24 
24 
24 


0.244(3) 
1.124(11) 
1.66(3) 


6.07(3) 
1.90(4) 
3.36(12) 


m'i 2.0 
m p 2.0 
mjv 2.0 


24 
24 
24 


0.100(6) 
0.94(2) 
1.38(2) 


6.31(6) 
2.10(6) 
3.65(8) 


m'i 2.0 


48 


0.047(10) 


6.47(11) 


dynamical with 


mfy^ = 0.06: 




m'i 5.325 


24 


0.332(3) 


5.16(2) 


m'i 2.0 


24 


0.120(8) 


6.44(8) 



3. DYNAMICAL QCD SPECTRUM 

To support our thermodynamics studies 0, 
we have measured the dynamical QCD hadron 
masses at zero temperature to set the scale. Un- 
less indicated otherwise, all the simulations dis- 
cussed in this section are done with lattice size 
8 3 x 32, m = 1.90, L s = 24, m [ * yn) = 0.02, 
and valence quark mass ranging from 0.02 to 0.22 
with an increment of 0.04. Valence extrapolations 
for some of the simulations discussed below are 
shown in Table |IJ 

Using DWF and the Wilson action, at Nt = 4 
the transition occurs at about (3 = 5.325. At 
this (3, we find m p = 1.18(3), m^ = 0.654(3) at 
m (d V n) = Q Q2j which givcg ^ = 163 ( 4 ) MeV an d 

m^/a = 427(ll)MeV. Using a larger 16 3 x 16 
volume only reduces m^ to 0.652(3), which sug- 
gests that the finite volume effect is very small. A 
Ward identity evaluation]^] shows that the resid- 
ual mass caused by the mixing between the two 
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walls plays an important role in this heavy pion 
mass. However, that study shows these residual 
mass effects do vanish in the limit of large L s for 
these full QCD simulations. 

We have also measured the masses for mj yn ^ — 
0.06 at j3 = 5.325. Performing real dynamical 
extrapolations using the two dynamical points 
obtained from to/ = 0.02,0.06, we get m% = 
0.320(6) + 5.38(ll)m^ n) . Compared with va- 
lence extrapolation, the dynamical extrapolation 
slightly decreases the pion mass as mj 0. 

As in the quenched studies, we have also in- 
vestigated the renormalization group improved 
gauge action. We have simulated at (3 = 1.90, 
ci = —0.331 which is about the transition point 
for Nt = 4 with DWF and the Iwasaki gauge ac- 
tion!]. We obtain m p = 1.16(2), = 0.604(3) 
at mf yn) = 0.02, which gives T c = 166(3)MeV 
and m^/a = 400(7)MeV. Surprisingly this is 
about the same value as that obtained at f3 c us- 
ing the Wilson action. Therefore, although the 
Iwasaki action helps to reduce the pion mass in 
our quenched study, this is not true for the dy- 
namical case. 

We have also measured the hadron masses at 
(3 = 2.0 for both mf yn) = 0.02, 0.06. As with the 
Wilson gauge action, we can draw the same con- 
clusion that the dynamical extrapolation m~ = 
0.088(10) + 6.9(2)mf yn) gives a slightly smaller 
pion mass as to/ — > 0. 

To study the finite L s effect, we are currently 
doing a simulation at (3 = 2.0, ci = —0.331, and 
L s = 48. Figure || shows the valence extrapo- 
lations of the pion, rho, nucleon mass values for 
L s = 24 and the pion mass for L s = 48. For 
L s = 24, we have m 7r (m^" ) = 0.02) = 0.475(7). 

For Ls = 48, we have obtained rn w {rri^ vn ^ = 
0.02) = 0.420(10). This confirms that the mixing 
between the two walls is at least a major cause of 
the heavy pion mass at /3 C , and this effect can be 
reduced by increasing L s . 

4. CONCLUSIONS 

Although for the quenched QCD spectrum the 
Iwasaki action lowers the pion mass compared 
with that obtained using the Wilson action, the 
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Figure 2. The quantities to^, m p , toat plotted 
versus m^ al ^ at (3 = 2.0, c\ — —0.331, too = 1.90, 
m ( d v«) _ q Q2 orL g3 x 32 using the Iwasaki action. 
The solid lines are valence extrapolations for to^., 
m p and to at at L s — 24. The dashed line is the 
valence extrapolation for at L s = 48. 

non-zero value of m. K {mf — > 0) is still problem- 
atic. From our dynamical QCD spectrum study, 
we have obtained a heavy pion mass at f3 c for 
Nt — 4 for both the Wilson and Iwasaki gauge 
actions at L s = 24. This can be improved by in- 
creasing the size of L s . These calculations were 
performed on the QCDSP machines at Columbia 
and RIKEN/BNL. 

REFERENCES 

1. D. Kaplan, Phys. Lett. B288 (1992) 342. 

2. V. Furman and Y. Shamir, Nucl. Phys. B439 
(1995) 54. 

3. P. Chen, et. al. Nucl. Phys. B (Proc. Suppl.) 
73 (1999) 204. 

4. T. Blum and A. Soni, these proceedings. 

5. Y. Iwasaki, Nucl. Phys. B258 (1985) 141; 
Univ. of Tsukuba report UTHEP-118 (1983), 
unpublished. 

6. P. M. Vranas, these proceedings. 

7. G. Fleming, these proceedings. 



